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The quark mean field model, which describes the nucleon usiug the constituent quark model, 
is applied to investigate the properties of finite nuclei. The couplings of the scalar and vector 
mesons with quarks are made density dependent through direct coupling to the scalar field so as 
to reproduce the relativistic Brueckner-Hartree-Fock results of nuclear matter. The present model 
provides satisfactory results on the properties of spherical nuclei, and predicts an increasing size of 
the nucleon as well as a reduction of the nucleon mass in the nuclear environment. 
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I. INTRODUCTION 



One of the most interesting topics in nuclear physics is to study the variations of hadron properties in nuclear 
' medium. Experimentally, the EMC effect reveals the medium modification of the internal structure of the nucleon [0. 
The EMC effect motivated many theoretical works on the study of hadrons in terms of quarks and gluons. At present, 
T-H ■ we are still far away to describe nucleons and nuclei in terms of quarks and gluons using quantum chromodynamics 
K>" ' (QCD), which is believed to be the fundamental theory of strong interactions. Thus, it is desirable to build models 
which incorporate quark and gluon degrees of freedom and respect the established theories based on hadronic degrees 
3^ of freedom. Such models are necessarily crude from various view points since the study of the nuclear many-body 
■ systems on the fundamental level is intractable. 

' The quark-meson coupling (QMC) model proposed by Guichon ^ provides a simple and attractive framework to 
^—1 incorporate quark degrees of freedom in the study of nuclear many-body systems. The model describes nuclear matter 
as nonoverlapping MIT bags in which the quarks interact self-consistently with structureless scalar and vector mesons 
, in the mean field approximation. The QMC model has been subsequently extended with reasonable success to various 
> - problems of nuclear matter and finite nuclei . The quark mean field (QMF) model as named in Ref . took the 
,—1 ' constituent quark model for the nucleon instead of the MIT bag model, so that it could overcome the shortcoming of 
^ : the QMC model as explained in Ref. The QMF model has been applied to study the properties of finite nuclei, 
where the coupling constants were determined by the saturation properties of nuclear matter |^ . 

In this paper, we introduce density dependent couplings for the QMF model, so that we can reproduce the relativistic 
Brueckner-Hartree-Fock (RBHF) results of nuclear matter Q. The RBHF approach as a microscopic many-body 
theory based on hadronic degrees of freedom, provided very closely the nuclear matter saturation properties, where 
the coupling constants and the form factors of the one-boson exchange potential (OBEP) were fixed from the nucleon- 
I nucleon scattering and the deuteron properties. The RBHF approach can be considered as parameter-free many-body 
theory, since all parameters in RBHF approach were adjusted to the two-nucleon problem at first. Encouraged by the 
success of the RBHF approach in describing nuclear matter properties, there were many applications of the RBHF 
results on finite nuclei [10|-^. The relativistic mean field (RMF) theory with density dependent couplings was 
developed for finite nuclei, as well as reproducing the RBHF results of nuclear matter. The RMF theory with density 
dependent couplings provided very good results on the nuclear properties as the binding energies and the nuclear 
radii with including the rearrangement contributions pl]-p^. Both the RBHF and RMF approaches are based on 
hadronic degrees of freedom, where the nucleons and mesons are treated as structureless particles. Therefore, they 
are not applicable to study the quark effect in nuclei. 

Hence, it is very interesting to connect the QMF model to the RBHF results by density dependent couplings, 
so that we can investigate the medium modification of the nucleon structure and finite nuclei properties, as well 
as reproducing the quantitative results of RBHF theory for nuclear matter simultaneously. The consideration to 
introduce density dependent couplings of mesons with quarks is based on the medium modifications of the meson 
properties. There are many works discussing the variations of hadron properties jl^-IT^ . In the present model, the 
scalar mean field provides a strong attraction and as a consequence the constituent quark mass is reduced in nuclear 
medium. The mesons, which are actually composites of quarks, are influenced due to the presence of the scalar mean 
field. Therefore, the density dependent couplings modeled through direct coupling to the scalar field represents, in a 
way, self-interaction in the scalar field. 
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The paper is organized as follows. In Sec. II we describe the QMF model with density dependent couplings. In 
Sec. Ill we present the parameters used in this paper, and compare the nuclear matter results with RBHF results. 
We then apply this model to finite nuclei and discuss the variations of the nucleon properties inside nuclei in Sec. IV. 
Section V is devoted to the summary. 



II. QUARK MEAN FIELD MODEL WITH DENSITY DEPENDENT COUPLING 

In this section, we give a brief description of the QMF model with density dependent couplings. We follow the basic 
concept of the QMF model [Qjs), in which the couplings of the mesons with quarks are constants. In the present work, 
we introduce density dependent couplings for the scalar (a) and vector (lu) mesons with quarks, which are modeled 
as functions of the a mean field. We include also the isovector meson, p, whose coupling with a quark is assume to 
be a constant for simplicity. We note that the contribution from the p meson is very weak for stable nuclei compared 
with those from the a and uj mesons. In addition, the photon field should be included for finite system. 

In the QMF model, the nucleon in nuclear medium is described in terms of the constituent quark model, in which 
the quarks couple with the a, u, and p mesons. We take the relativistic mean field (RMF) approximation |I|-|2^ for 
the meson fields. In the constituent quark model, the quarks in a nucleon satisfy the following Dirac equation: 

[*7m9'' - nig + gl{a)a - g^^{a)uj^" - g^^prsj" - Xc] q{r) = 0. (1) 

where niq is the constituent quark mass and Xc the confinement. The g^ic) and g^icr) as functions of the a mean 
field, are the couplings of the a and lo mesons with quarks, while the g'^ denotes the coupling constant of the p meson 
with a quarks. We follow Ref. to take into account the spin correlations and remove the spurious center of mass 
motion, and then work out the effective nucleon mass M*. 

For a static nuclear many-body system, the Lagrangian within the mean field approximation can be written as 

C = ^ [tj^d^^ - M*{a) - gU^Pl"" - 5pPT37°] ^ (2) 

-1 [(VCT)^ + mla'] + \ [(V-)^ + mlu'] + \ [{VP? + my] 
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Here, ■0 denotes the nucleon fields; cr, and p are the cr, w, and p meson mean fields with masses mo-, m^, and 
TTip, respectively; A denotes the photon field. The nuclear many-body system can be solved with the change of the 
nucleon properties under the influence of the scalar mean field, which is exclusively expressed in the effective nucleon 
mass, M*(g), as a function of the a mean field. The uj and p mean fields do not cause any change of the nucleon 
properties pl[ , and they appear merely as the energy shift. The relations between nucleon- meson couplings and 
quark- meson couplings have been given in Ref. j|] as gu){cr) — 3g^{(T) and gp — gf,. We will apply the present model 
to study the properties of nuclear matter and finite nuclei in the following two sections. 



III. PARAMETERS AND NUCLEAR MATTER RESULTS 



In this section we will show the numerical results of nuclear matter and discuss the density dependent couplings. 
The change of the nucleon properties under the influence of the a mean field has been studied in Ref. ||], where the 
effective nucleon mass, M*, has been shown as a function of the quark mass correction, SrUq. Here Sruq is related 
to a mean field as Sruq = g^{a)a. We take the same quark mass and two types of confinement (xc ~ ^kr'^ and 
Xc = ifcr2(l + 70)/2 with k = 1000 MeV/fm^) as those used in Ref. @. 

To perform the nuclear matter calculation, we take the Lagrangian given in Eq. (2), in which all the derivative 
terms of the meson fields vanish for infinite nuclear matter due to the translational invariance. Then, the nucleons 
and mesons obey the following Euler-Lagrange equations, 

[tj^d^ - M:(a)-.9^(aV7°-5pPT37°]^ = 0, (3) 



ml a 



dM*{a) 



da 

dM*{6mq) 
dSmo 



dgu{<y) 



da 



'(^7V) 



dM:{5mq) dg%{a) 



dSmo 



da 



(4) 
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^ Lo = .g^(a)(^7V) =35^(^)(Vi7», (5) 



"^-p P = gp(?AT-37''V') = 3p(V'T37°'0)- (6) 



Here, the bracket ( ) means the expectation value of the operator between the nuclear ground state. We note that 
the appearance of the last two terms in Eq. (4) is due to the variable couplings, which is similar to the inclusion of 
the rearrangement terms in the density dependent RMF approach based on hadronic degrees of freedom Q. The 
effective nucleon mass M* and its derivative with respect to the Sruq has been studied in Ref. |^ . We take the meson 
masses as = 550 MeV, = 783 MeV, and rup = 770 MeV. 

We introduce density dependent couplings, g%{(j) and g^{cr), in the present model. This consideration is based on 
the variation of meson properities as the nuclear environment changes. The mesons, which are actually composites of 
quarks, are influenced by the a mean field |p^ , |l7| . The quark- meson couplings in nuclear medium may be modified as 
density increases. It is therefore reasonable to attempt density dependent couplings, which contain the complicated 
medium modification effectively. This treatment is similar in spirit to the density dependent RMF approach based 
on hadronic degrees of freedom, which provided a realistic description of bulk properties of finite nuclei and nuclear 
matter [p^|-p^. The nucleon-meson couplings in the density dependent RMF approach were found to decrease with 
increasing nuclear matter density. This suggests a decreasing quark-meson couplings due to the relations between the 
nucleon-meson couplings with the quark-meson couplings. 

we model the density dependent couplings, and g^, through direct couplings to the a mean field. We parameterize 
these variable couplings as 

glia) =a.~ b,a'l\ (7) 

gl{^) = ci^-K<y^"'- (8) 

With these formulas, it is possible to reproduce the scalar and vector potentials of the RBHF theory for symmetric 
nuclear matter within a good accuracy. We list the parameters of and hi in Table I, which are determined by the 
RBHF results with potential A Q. The parameter 5^, which is still a constant for simplicity, is taken as gj, = 4, so 
that we can get the symmetry energy to be around 35 MeV. 

In Fig. 1, we plot the scalar and vector potentials. Us and f/y, as functions of the nuclear matter density p. The 
RBHF results Q are marked by solid dots, as well as the results in the RMF model with TMI parameter set by 
dotted curve for comparison. We show in Fig. 2 the energy per nucleon, Ej as a function of density p. The solid 
and dashed curves, which are the results of the present model with two types of confinement, are almost identical. 
We note that the saturation density, binding energy and incompressibility of the solid curve in Fig. 2 are 0.16 fm~^, 
— 15.5 MeV and 159 MeV, respectively. We found there exist slight differences between the RBHF results and those 
in the present model in Fig. 2, which is partly due to the appearance of the rearrangement terms in Eq. (4). The 
importance of rearrangement terms in a relativistic density dependent field theory has been discussed extensively in 
Ref. [O. 



IV. FINITE NUCLEI RESULTS 



In this section we will present the calculated results of some spherical nuclei and discuss the variations of the 
nucleon properties inside nuclei. The meson mean fields in a spherically symmetric nucleus are functions of the radial 
coordinate r. Here, we follow the prescription of Ref. [^,^ and assume the nucleons obey the Dirac equation. 



i-ipd^ - Ml{o) - g^{o)u:{r)i' - gpp{r)T^^'' - -(1 + T^)A{r)i' 



(9) 



where the effective nucleon mass is a function of the scalar mean field g . The meson mean fields (cr, u^, p) and 
the photon field (A) satisfy the following Klein-Gordon equations, 



dSm„ 



dSnin 



da 



a(r)(^^)+3^^c.(r)(^7V), 



(10) 
(11) 
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Ap(r) - w?pp(r) = 



(12) 



A^(r) = -^(V;(l + r3)7V). (13) 

We solve the above equations self-consistently using the parameters determined by the RBHF results of nuclear 
matter with potential A. In Fig. 3 and 4, we plot the charge density distributions for '^''Ca and ^"^Pb, respectively. 
The experimental values taken from Ref. are shown by solid curves. The calculated results with scalar confining 
potential (xc = ^kr'^) are shown by dashed curves, while those with scalar-vector confining potential [xc — ^kr^{l + 
7^)/2) by dash-dotted curves. The curves with two types of confinement are very close, since the parameters used here 
are determined by the same RBHF results. The theoretical central densities are slightly higher than the experimental 
values, which is related to the RBHF results used in determining the parameters. The dependence of the finite nuclei 
properties on the RBHF results were discussed in more detail in Refs. ||l^,|l^. The calculated results in the present 
model for the binding energies per nucleon E/A^ the rms charge radii i?c, and the spin-orbit spitting /\Els are 
compared with the experimental values in Table H. The experimental data taken from Refs. ||2^,Q are listed in the 
last column. The results with two types of confinement agree very well with the experiment, and no explicit difference 
between them. In Table HI and IV, the single particle energies of ''"Ca and ^"®Pb are compared with the experimental 
data IUJ. We show in Fig. 5 and 6 the scalar and vector mean fields {a and uj) as functions of r for ''"Ca and ^"^Pb. 
Here, the results with Xc = \kr^ are given by solid curves, while those with Xc = \kr'^{l + 7°)/2 by dashed curves. 
Again, the corresponding results with two types of confinement are very similar due to the same reason as mentioned 
before. 

It is very interesting to see the variations of the nucleon properties inside nuclei. We plot in Fig. 7 and 8 the 
ratios of the nucleon rms radius and effective mass in ^"Ca and ^osp]^ those in free space, R/Rq and M*/AIn, as 
functions of nuclear radial coordinate r. We find the model with Xc = ^kr^{l + 'j'^)/2 (dashed curves) gives near 
10% increasing nucleon radius at the center of the nuclei, while that with Xc = ^kr^ gives around 5% increase. The 
nucleon radius decrease to its value in free space from center to outside, as the scalar mean field decreasing. As for 
the effective nucleon mass, both scalar and scalar-vector confinement cases give around 40% reduction at the center 
of the nuclei. The curves with two types of confinement are almost identical. It is due to the reproduction of the 
RBHF scalar potential. 

V. CONCLUSION 

We have developed the QMF model with density dependent couplings, so that we can investigate the variations 
of the nucleon structure and finite nuclei properties, as well as reproducing the RBHF results of nuclear matter 
simultaneously. We introduced the density dependence of the couplings of mesons with quarks through direct couplings 
to the (T field based on the capability of a field in influencing hadron properties. We have used the constituent quark 
model for the nucleon, which naturally allows the direct coupling of cr, and p mesons with quarks. With the 
parameters determined by the RBHF results of nuclear matter, we can perform the numerical calculations for finite 
nuclei. 

We have found the QMF model with density dependent couplings gave very successful descriptions of spherical 
nuclei. By comparison with the density dependent RMF approach on hadronic level, we have reproduced the properties 
of finite nuclei with the same or even better agreement to the experimental data. The present model with scalar 
confinement has the capability to produced better results on the binding energies per nucleon compared with those in 
Ref. Very weak dependence on the type of the confinement was found in the results of nuclear properties, which 
is due to the reproduction of the RBHF scalar and vector potentials for nuclear matter. 

We have also investigated the variations of the nucleon properties in spherical nuclei. The QMF model with density 
dependent couplings can provide a significantly swollen nucleon radius in nuclear environment. We found the model 
with Xc — ^kr^{l + 7°)/2 predicted near 10% increasing nucleon radius at the center of nuclei, while that with 
Xc = ^kr^ gave around 5% increase. At the same time, around 40% reduction of the nucleon mass at the center of 
nuclei were predicted. 

ACKNOWLEDGMENTS 

One of the authors (H.S.) would like to thank H. Toki for valuable discussions. This work was supported in part 
by the National Natural Science Foundation of China. 



4 



[1] J. J. Aubert et al, EMC Collaboration, Phys. Lett. 123B, 275 (1983). 

[2] P. A. M. Guichon, Phys. Lett. 200B, 235 (1988). 

[3] K. Saito and A. W. Thomas, Nucl. Phys. A574, 659 (1994). 

[4] P. A. M. Guichon, K. Saito, E. Rodionov, and A. W. Thomas, Nucl. Phys. A601 349 (1996). 

[5] H. MuUer and B. K. Jennings, Nucl. Phys. A640, 55 (1998). 

[6] P. Wang, R. K. Su, H. Q. Song, and L. L. Zhang, Nucl. Phys. A653, 166 (1999). 

[7] H. Toki, U. Meyer, A. Faessler, and R. Brockmann, Phys. Rev. C 58, 3749 (1998). 

[8] H. Shen and H. Toki, Phys. Rev. C 61, 045205 (2000). 

[9] R. Brockmann and R. Machleidt, Phys. Rev. C 42, 1965 (1990). 

[10] R. Brockmann and H. Toki, Phys. Rev. Lett. 68, 3408 (1992). 

[11] C. Fuchs, H. Lonske, and H. H. Woltcr, Phys. Rev. C 52, 3043 (1995). 

[12] F. Ineichen, M. K. Weigel, and D. Von-Eiff, Phys. Rev. C 53, 2158 (1996). 

[13] H. Shen, Y. Sugahara, and H. Toki, Phys. Rev. C 55, 1211 (1997). 

[14] G. E. Brown and M. Rho, Phys. Rev. Lett 66, 2720 (1991). 

[15] C. Song, G. E. Brown, D. -P. Min, and M. Rho, Phys. Rev. C 56, 2244 (1997). 

[16] K. Saito and A. W. Thomas, Phys. Rev. C 51, 2757 (1995). 

[17] K. Saito, K. Tsushima, and A. W. Thomas, Phys. Rev. C 55, 2637 (1997). 

[18] B. D. Scrot and J. D. Walccka, Adv. Nucl. Phys. 16, 1 (1986). 

[19] Y. K. Gambhir, P. Ring, and A. Thimet, Ann. Phys. (N.Y.) 198, 132 (1990). 

[20] Y. Sugahara and H. Toki, Nucl. Phys. A579, 557 (1994). 

[21] S. Fleck, W. Bcntz, K. Shimizu, and K. Yazaki, Nucl. Phys. A510, 731 (1990). 

[22] Y. Shibata and H. Tezuka, Z. Phys. C 62, 533 (1994). 

[23] H. de Vries, C. W. de Jager, and C. de Vries, At. Data Nucl. Data Tables 36, 495 (1987). 

[24] X. Campi and D. W. Sprung, Nucl. Phys. A194, 401 (1972). 



5 



TABLE I. Parameters needed to reproduce the RBHF results of nuclear matter with potential A in Ref. 





QMF 


QMF 






(Xc = ifcr2(l+7")/2) 




6.415 


7.328 




5.852 


6.016 




7.543 


7.543 




6.470 


6.470 



TABLE II. Binding energies per nucleon (E/A), rms charge radii (Re), and spin-orbit splittings for neutrons {AE^ 
protons (AEp) are compared with the experimental data [ESlpd]. The energies are in MeV, and the radii are in fm. 



and 







QMF 

iXc = \kr^) 


QMF 

(Xc = ifcr2(l + 7°)/2) 


Expt. 


^"Ca 


E/A 


8.65 


8.62 


8.55 




Rc 


3.46 


3.45 


3.45 




A£;„(ld5/2-ld3/2) 


4.5 


4.5 


6.3 




AEp{ld^/2 - ld3/2) 


4.4 


4.5 


7.2 


*«Ca 


E/A 


8.58 


8.55 


8.67 




Rc 


3.49 


3.48 


3.45 




AS„(ld5/2 - ld3/2) 


3.8 


3.8 


3.6 




AEj,[U5/2 - 1^3/2) 


3.9 


3.9 


4.3 




E/A 


8.50 


8.49 


8.71 




Rc 


4.29 


4.28 


4.26 




AS„(2p3/2 - 2pi/2) 


1.2 


1.3 


0.5 


208 p^^ 


AEp(2p3/2 - 2pi/2) 


1.2 


1.2 




E/A 


7.62 


7.62 


7.87 




Rc 


5.57 


5.57 


5.50 




Ai5„(2/7/2- 2/5/2) 


1.6 


1.6 


1.8 




^Ep{lgs/2 - 157/2) 


2.6 


2.5 


4.0 



TABLE III. Single particle 
All energies are in MeV. 


energies 


of neutron (n) 


and proton (p) for ^"Ca. The experimental data 


are taken from Ref. . 


SheU 


n 


QMF 

= \kr^) 
P 


(xc = i 

n 


QMF 

fcr2(l + 7»)/2) 
P 


Expt. 
n 


P 


lSl/2 
lP3/2 
lPl/2 
1^5/2 
1^3/2 
2Sl/2 


51.3 
36.2 
33.6 
21.4 
16.9 
17.3 


43.1 
28.4 
25.7 
14.0 
9.5 
10.0 


51.5 
36.3 
33.7 
21.5 
17.0 
17.3 


43.3 
28.5 
25.8 
14.0 
9.6 
10.0 


50.0 
30.0 
27.0 
21.9 
15.6 
18.2 


50±11 
34±6 
34±6 
15.5 
8.3 
10.9 
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TABLE IV. Single particle energies of neutron (n) and proton (p) for ■^''**Pb. The experimental data are taken from Ref. 
All energies are in MeV. 



QMF QMF Expt. 



Shell 


n 


P 


n 


P 


n p 


lSl/2 


56.3 


43.4 


56.1 


43.4 




lP3/2 


50.2 


38.0 


50.1 


37.9 




lPl/2 


49.9 


37.6 


49.8 


37.5 




1^5/2 


43.1 


31.3 


43.0 


31.3 




1^3/2 


42.3 


30.4 


42.2 


30.3 




2Sl/2 


39.7 


27.3 


39.7 


27.3 




1 /7/2 


35.2 


23.7 


35.1 


23.7 




1/5/2 


33.6 


22.1 


33.6 


22.1 




2P3/2 


30.2 


18.2 


30.2 


18.2 




2Pl/2 


29.6 


17.5 


29.6 


17.5 




I59/2 


26.6 


15.5 


26.6 


15.5 


15.4 


I57/2 


24.2 


13.0 


24.2 


13.0 


11.4 


2^5/2 


20.9 


9.1 


20.9 


9.2 


9.7 


2d3/2 


19.7 


7.9 


19.7 


8.0 


8.4 


3Sl/2 


18.9 


6.8 


18.9 


6.9 


8.0 


lftll/2 


17.8 


6.9 


17.8 


7.0 


9.4 


1^9/2 


14.4 




14.4 




10.8 


2/7/2 


11.9 




11.9 




9.7 


2/5/2 


10.3 




10.3 




7.9 


3p3/2 


9.6 




9.6 




8.3 


3pi/2 


9.0 




8.9 




7.4 


1*13/2 


8.9 




8.9 




9.0 
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FIGURE CAPTIONS 



Figure 1: The scalar and vector potentials of nuclear matter, Us and C/y, as functions of the nuclear matter density p. 
The results of the present model with Xc — ^kr^ are shown by solid curves, while those with Xc = ^A:r^(H-7°)/2 
are shown by dashed curves. The RBHF results with potential A in Ref. ||] are marked by solid dots, and the 
results in the RMF(TMl) model are plotted by dotted curves for comparison. 

Figure 2: The energy per nucleon, E/A, as functions of the nuclear matter density p. The curves are labeled as in 
Fig. 1. 

Figure 3: The charge density distributions for ^"Ca compared with the experimental data (solid curve) The 
dash-dotted and dashed curves are the results in the present model with Xc — \kr'^ and Xc = ^kr^ll + 7*')/2, 
respectively. 

Figure 4: Same as Fig. 3 but for ^"^^Pb. 

Figure 5: The scalar and vector mean fields, a and cj, as functions of the radial coordinate r for ^"Ca. The results in 
the present model with Xc = \kr^ are shown by solid curves, while those with Xc — ^kr^{l + j'^)/2 are shown 
by dashed curves. 

Figure 6: Same as Fig. 5 but for ^"^Pb. 

Figure 7: The ratios of the nucleon rms radius and effective mass in ^"Ca to those in free space, R/Rq and M*/Mn, 
as functions of the radial coordinate r. The results in the present model with Xc = ^kr"^ are shown by solid 
curves, while those with Xc = jkr^{\ + 7°)/2 are shown by dashed curves. 

Figure 8: Same as Fig. 7 but for ^''^Pb. 
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